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Abstract 

The goal of this work is to investigate the perturbation of the climate system due to 
large impact events. 

Impacts are among the most important mechanisms for the evolution, distribution, and 
destruction of life in the universe. However, the possible climatic effects of an impact 
were not seriously considered until 1980, when Louis and Walter Alvarez suggested that 
the profound end-Cretaceous extinction might have been caused by the impact of an 
asteroid or comet about 10 km in diameter. Since then, the climatic change associated 
with the end-Cretaceous impact has become one of the most interesting and still 
unresolved questions in linking the well-known Chicxulub impact event and the end- 
Cretaceous mass extinction. While the end-Cretaceous impact offers the best-documented 
case of an impact affecting the Earth’s climate and biota, even smaller (and more 
frequent in time) impacts could introduce significant perturbations of the climate 
comparable, if not larger, to the largest known volcanic perturbations. 

We propose to study the mechanical and thermal state of the atmosphere following an 
impact event. This will be done by using both one-dimensional and three-dimensional 
climate models. When necessary, modifications of the state-of-the-art general circulation 
models will b e carried out. 

We want to use the end-Cretaceous impact event as a case study. This allows us to 
take advantage of the extensive modeling of this impact event that has already been 
carried out through a previous Exobiology grant (NAGW-5159). Furthermore, a large 
experimental dataset, that can be used to constrain and test our models, is associated with 
the end-Cretaceous mass extinction (one of the largest of the Phanerozoic) and impact 
event. 

1. Research Results 

In this project we have focused on characterizing climate forcing (change in the 
Earth’s total radiative balance at the tropopause; [W/m 2 ]) and climate sensitivity (mean 
change in global surface temperature occurring in response to a specific radiative forcing; 
[K.]) of impact-related loads of S-bearing gases into the stratosphere using one- 
dimensional (ID) atmospheric models from the National Center for Atmospheric 
Research (NCAR, Boulder, CO) coupled to a simple sulfate aerosols model. Some of the 
results have been summarized in a paper that was published in Astrobiology as part of the 
Rubey colloquium proceedings (Pierazzo et al., 2003). This paper also serves as a general 
review of modeling work done in the past to understand abrupt climate perturbation 
associated with large impact events and nuclear wars. 

Results of forcing and sensitivity studies using l D atmospheric models include; 
Saturation Effect : Estimated K/T impact-related injection of S-bearing gases from the 
sedimentary layer range from about 20 to 100 Gt, causing a negative forcing of a few 
hundred W/m 2 , about two orders of magnitude stronger than the Pinatubo volcanic 



eruption, for up to two years after Fig. 1: Climate forcing (cooling) for various 
the impact. The simulations indicate 
a " saturation effect", Fig.l: the 
magnitude and duration of the 
forcing does not depend on the 
amount of S injected for loads 
above about 25-30 Gt. The forcing 
due to S-loading from the projectile 
alone is around one order of 
magnitude larger than for Pinatubo. 

Rough estimates of the forcing due 
to injection of dust (Fig. 1 ; Covey et 
al., 1994) are about one order of 
magnitude inferior to those relative 
to S-loading. For comparison, we 
estimated that the positive forcing 
(heating) from the impact-related 
injection of CO: is -1-3 W/m 2 , 
similar to the estimated forcing 
from greenhouse gases due to 
industrialization. 

Residence Time : The stratospheric 
residence time of gases and 
particulate affects the duration of 
the forcing. Our initial assessment of climate forcing associated with the impact event 
assumed a stratospheric residence time typical of a steady-state atmosphere that is 
around 2 years. However, in an impact-perturbed atmosphere heating of the upper 
atmosphere and cooling of the surface may cause a stronger atmospheric stratification 
and favor smaller aerosol sizes. In this case, the stratospheric residence time for tracer 
gases ad sulfate aerosol may be much larger than in the unperturbed atmosphere. This 
increases the overall duration of the impact- related effects. Tests with the Sulfate 
Aerosol Model (SAM) coupled to the Column Radiation Model (CRM) from NCAR 
indicate that the main effect of increasing the stratospheric residence time by a factor 
of 2 to 5 (corresponding to between 4 and 10 years) from the adopted baseline value 
is a significant increase in the duration of the impact-related effects but hardly affects 
the maximum change in the net radiative fluxes at the tropopause (i.e., the radiative 
forcing). 

Aerosols’ Characteristics : In this study we use a constant sulfate aerosol density of 1.83 
g/cm 3 (corresponding to a pure sulfate aerosol). In reality, the density of sulfate 
aerosols depends on the sulfate-water solution, which, in turn, varies with 
temperature. We found that varying the density between 1.5 g/cmg 3 (equivalent to a 
solution of 60% H 2 SO 4 and 40% H 2 0) and 2 g/cm 3 (due, for example, to the addition 
of impurities to the sulfate aerosols) introduces at most an uncertainty of 15% in the 
forcing. The aerosols optical parameters (extinction coefficient, single scattering 
albedo, asymmetry parameter, and forward scattering fraction) needed for the 
atmospheric models are estimated using Mie theory. Small impurities, especially dust 


stratospheric S-loadings (20% S0 3 , 80% SO;). 
Solid lines: 10, 30, 100, 300 Gt; dashed lines (2 
Gt): Maximum loading from the projectile. 
Modeled forcing form a Pinatubo-like injection 
(20 Mt SO;, dot-dashed line) and observed 
Pinatubo forcing ( diamonds ; McCormick et al., 
1995) are shown for comparison. Gray line : dust 
loading (Covey et al., 1994). 
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and soot, strongly affect the optical properties of the aerosols, especially in the visible 
and near-infrared, by modifying the imaginary part of the index of refraction. This in 
turn affects the transmission in the aerosol layer. The overall effect, however, does 
not appear to be dramatic, and it has few repercussions on the overall radiative 
forcing. 

S-bearing Gases : While there is little doubt that significant amounts of S-bearing gases 
were released in the impact event, the distribution of S among various S-oxides is not 
well characterized (e.g., Ohno et al.. 2003). Our results indicate that the release of S 
mainly as SO 3 causes an early production of sulfate aerosols, thus maximizing the 
abruptness and magnitude of the climate forcing. The release of sulfur as SO 2 slows 
down the sulfate-formation process, resulting in a smaller initial forcing that extends 
overall over a longer period. In particular we find that although the overall aerosol 
production for the SO 2 case is almost two orders of magnitude smaller than that of the 
SO 3 case the maximum radiative forcing decreases to only about two-thirds of the 
SO 3 case (in agreement with the saturation effect discussed in point a) ). 

Climate Sensitivity : As a next step in complexity, we coupled the SAM to NCAR's 
single column model, SCCM. Although this does not provide a complete assessment 
of the mean change in global surface temperature in response to the atmospheric S- 
loading, it allows us to investigate the immediate response of the atmosphere to 
perturbation from impact-related sulfate aerosols. SCCM is equivalent to a grid 
column of the more complete global climate model CCM3 (Kiehl et al., 1996) where 
the performance of the parameterized physics for the column is evaluated in isolation 
from the rest of the large-scale model. While lacking the more complete feedback 
mechanisms available to an atmospheric column imbedded in the global model, it 
provides an inexpensive first look at the response of the system to the forcing 
introduced by a particular parameterization. The lack of horizontal feedbacks, 
however, implies that we can only evaluate climate sensitivity immediately following 

the injection of the S-bearing Fjg. 2 : Combined impact-related S (100 Gt) and CO : 
gases, up to about a month or two 
after the impact loading. Any 
longer-term climate sensitivity 
can only be investigated with a 
fully coupled climate system 
model. 

The presence of the S-bearing 
gases and sulfate aerosols (strong 
LW absorbers) initially in the 
upper atmospheric layer of the 
model produces a significant 
change in the atmospheric 
radiation fluxes. This results in a 
strong heating of the stratosphere 
accompanied by a strong cooling 
at the Earth's surface (Pierazzo & 

Hahmann, 2001; Pierazzo et al., 

2003). Compared to a Pinatubo- 


loading at 37°N, 97°W (North America) for present 
(pre-industrial= 280 ppm; solid line) and end- 
Cretaceous atmospheric CO, (2-' pre-industrial- 580 
ppm; 4><pre-industrial= 1 120ppm). 
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type eruption, the model estimates that in the uppermost layer the temperature 
increases by (at least) several tens of degrees, more than an order of magnitude that 
associated with Pinatubo. At the surface the impact-produced cooling is around 
several degrees, again at least an order of magnitude that associated with Pinatubo. 
Temperature changes are strongest over continental interiors than over the oceans, 
regardless of latitude or season. Even at coastal location the nearby vicinity of the 
oceans act to mitigate the temperature changes due to the stratospheric aerosols. At 
the top of the atmosphere, however, the influence of the oceans below is minimal, and 
latitudinal (high versus low) and seasonal (summer versus winter) effects become the 
main factors influencing the magnitude of temperature changes due to stratospheric 
aerosols. Adding the impact-loading of CO 2 to the S-loading effect mitigates only 

slightly the overall effect, as 
Fig. 3: Zonal average temperature deviations from the shown in Fig. 2. 
unperturbed case for a 0.1, 2, and 10 Gt S-loading in 

CCM3^est runs° SPhere 8 ^ ***' injeCti0 "’ *° m CCM3 testing : We have begun 
* testing of NCAR’s three- 



dimensional (3D) atmospheric 
model CCM3 (Kiehl et al., 1996) 
after coupling the sulfate aerosol 
model. Contrarily to the ID 
atmospheric codes, the 3D code 
has continuously been updated to 
improve the various physics 
components and code coupling. 
Therefore, the coupling of the 
SAM to the 3D code is not the 
same as for the 1 D code, and had 
to be further tested and refined. 
Preliminary tests have highlighted 
few problems that have helped in 
improving the overall code. 
Figure 3 shows CCM3 test runs 
results for various S-loads before 
the 10 Gt case blows up (8days). 
The size of perturbations 
associated with the impact is 
already revealing limitations of 
some parameterizations used in 
the atmospheric model. For 
example, in CCM3 the 
parameterization of the saturation 
vapor pressure of water is valid 
for temperatures up to 375. 16K; 
however, impact-effects can cause 
stratospheric temperature 

increases above that maximum 
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limit, causing the model to crash. This parameterization must be adjusted to allow for 
a wider temperature range. It is also becoming clear that the vertical resolution of 
CCM3 in the stratosphere is too coarse. This makes the model less stable when very 
large perturbations are involved. A higher vertical resolution in the stratosphere may 
help stabilize the simulation all together, but this requires a different atmospheric 
model and the availability of large and fast computers. 
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